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Simple models of earthquake faults are important for understanding the mechanisms for their 
observed behavior, such as Gutenberg-Richter scaling and the relation between large and small 
events, which is the basis for various forecasting methods. Although cellular automaton models 
have been studied extensively in the long-range stress transfer limit, this limit has not been studied 
for the Burridge-Knopoff model, which includes more realistic friction forces and inertia. We find 
that the latter model with long-range stress transfer exhibits qualitatively different behavior than 
both the long-range cellular automaton models and the usual Burridge-Knopoff model with nearest 
neighbor springs, depending on the nature of the velocity-weakening friction force. This results have 
important implications for our understanding of earthquakes and other driven dissipative systems. 

PACS numbers: 91.30.Px, 02.60.Cb,05.20.-y, 05.45.-a 



Earthquake faults are examples of driven dissipative 
systems Models of fault systems are important for 
understanding scaling laws, the occurrence of character- 
istic events, and the relation between small and large 
earthquakes 0, Q, 0]- In addition to the benefits that 
would result from understanding earthquake faults, un- 
derstanding driven dissipative systems is important in 
physics and related areas. 

Simulations of the Burridge-Knopoff model .4] have led 
to much insight. The model consists of blocks connected 
by linear springs to their nearest neighbors with spring 
constant kc- The blocks are also connected to a loader 
plate by linear springs with spring constant k^, and rest 
on a surface with a nonlinear velocity-weakening stick- 
slip friction force that depends on a parameter a which 
controls how quickly the friction force decreases as the 
velocity is increased and determines the amount of stress 
dissipated in an event. The model was studied numeri- 
cally in one dimension Q and more recently S • 

An earthquake event is defined as a cluster of blocks 
that move due to the initial slip of a single block. The 
moment of an event is oc Uj , where the sum is over 
all the blocks in an event and Uj is the net displacement 
of block j . The main result of prior studies 0, 0, 01 is 
that the moment distribution, P{M), scales as for 
small localized events with an exponent b ~ 2 for a > 1. 
Large events show a pronounced peak in P{M) indicat- 
ing a significant presence of characteristic or system- wide 
events. In contrast, the number of earthquake events 
with s blocks, P{s), does not show power law scaling. 
Note that the blocks are connected only to their near- 
est neighbors in the Burridge-Knopoff model studied in 
Refs. IIIEIEH 

In the cellular automaton (CA) versions of the 



Burridge-Knopoff model H, with nearest-neighbor 
stress transfer, P{s) does not exhibit power law scal- 
ing A generalization of the CA models with 
long-range stress transfer yields considerable differences 
with the nearest-neighbor CA models |^| and with the 
nearest-neighbor Burridge-Knopoff model. In particu- 
lar, the long-range CA models run in a Gutenburg- 
Richter scaling mode and exhibit events of all sizes, con- 
sistent with the system being near a mean-field critical 
point Ila . [ij) ■ Quasi-periodic system- wide events 
are not observed, and P{s) exhibits power law scaling. 
Small and medium size events can be interpreted as fluc- 
tuations about a free energy minimum Large 
events drive the system out of "equilibrium" from which 
the system decays back to an equilibrium state . We 
refer to this behavior as punctuated ergodicity, which 
also has been observed in the Southern California fault 
system |T^ . 

Both the CA and nearest-neighbor Burridge-Knopoff 
models lack several elements that would make them more 
realistic representations of earthquake systems. In partic- 
ular, the latter does not include long-range stress trans- 
fer, and the long-range CA models do not have inertia 
and more realistic friction laws. 

In this letter we discuss the results of adding more re- 
alistic long-range stress transfer to the Burridge-Knopoff 
model. Our primary result is that the statistical proper- 
ties of this generalized Burridge-Knopoff model depend 
not only on the range R of the stress transfer, but also on 
the characteristics of the friction force. For large velocity 
weakening, that is, large a (see Eq. 0), the properties 
of the Burridge-Knopoff model are largely independent of 
R, although there are some characteristics of the model 
that differ from the nearest-neighbor {R—1) model and 
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are reminiscent of the long-range CA models. For small 
a and large i?, the results are identical to the long-range 
CA models. 

The usual {R — 1) Burridge-Knopoff model in one di- 
mension is governed by the equation of motion 



2xj + Xj^i) - kLXj - F{v + ij), (1) 



where Xj is the displacement of block j from its equilib- 
rium position, V is the speed of the substrate moving to 
the left with a fixed loader plate, F{x) = Fq4){x/v) is 
the velocity-dependent friction force, w is a characteristic 
velocity, and m is the mass of a block. 

As in Ref. 0] we introduce scaled variables r = ujpt^ 
— kL/m, and Uj = {kL/ FQ)xj, and rewrite Eq. (^) as 



ilj = ^^(uj+i — 2uj + Uj-i) — Uj — (f){2av + 2auj)^ (2) 

with 2a = ajpFo/kLV, P — k^/k^, and v — vk^ / {ujpFo)] 
the dots now denote differentiation with respect to r. 
The friction force is 



(-00,1], y^Q 
y>{). 
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We generalize the Burridge-Knopoff model by assum- 
ing that a block is connected to R neighbors (in each 
direction) with spring constant kc/R] R = 1 corresponds 
to the usual Burridge-Knopoff model. We use the second- 
and fourth-order Runge-Kutta algorithms with time step 
At = 0.001 to solve Eq. modified for arbitrary R, 
with the friction force given in Eq. (Q. Both algorithms 
and other fourth-order algorithms give similar results. 

A block is defined to be stuck if its speed is less than 
the parameter vq, its velocity is decreasing, and the stress 
(the force due to the springs coupled to it) is smaller than 
the maximum static friction force Fq (taken to be unity) . 
If a block is stuck, its velocity is set equal to zero, and 
the static friction force cancels the stress. A stuck block 
will accelerate if the stress is greater than Fq ; in this case 
the friction force drops to the value 1 — cr. An earthquake 
begins with the slip of a block and ends when all blocks 
become stuck. Blocks can become stuck and then slip 
again during an event. The results in this paper are for 
£ = 10, cr = 0.01, vo = 10-^ TV = 5000 blocks, i? = 1, 
100, 500, 1000, a = 0.0, 0.5, and 2.5, and 10^ events. 

We initially set iij = for all j and assign random 



displacements to all the blocks, 
on all the blocks and update iij 



We compute the force 
and Uj for all j. 



We 

continue these updates until all blocks become stuck. We 
then move the substrate (the loader plate is fixed) until 
the stress on one block exceeds unity. This stress loading 
mechanism is known as the "zero velocity limit" ^ and 
is equivalent to setting 1^ = 0. 

In Fig. n we plot P{M) for various values of R and a. 
For a = 2.5 and i? = 1, P(M) ~ M'^ with 6 « 2 as 
in Ref. 0. For larger values of R, b remains w 2, but 
the scaling range decreases and the characteristic earth- 
quakes (M > 1) become better defined. For a = 0.5 and 



R = 100, b = 2.01 
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FIG. 1; Distribution of the moment M of the events P(M). 
For Q = 2.5 and R = 1 the slope & « 2, consistent with 
Ref. 5]; 6 does not change significantly as R increases. For 
a — 0.5 the scaling range becomes better defined as R in- 
creases and the slope converges to b ~ 1.5. For a = 0, fe ~ 1.5 
for R > 100; there is no power law behavior for R — 1 and 
small Q. All plots have been displaced vertically for clarity. 

R > 100, b appears to approach 1.5 as R increase s, and 
there are fewer characteristic earthquakes (see Fig. 1(b) 1. 
The results in Fig. |l(c)| for a = are consistent with the 
long-range CA models jH m for which 6 = 3/2. The 
distribution of events, P{s), does not exhibit power law 
behavior for i? = 1 and a = 2.5, consistent with Ref. |^, 
but does so for larger R with P{s) ~ and 6 « 2 (see 
Fig. ID). For small a, b ^ 1.5 for R > 100. 

We now consider the existence of ergodicity and deter- 
mine the metric fl / (t) ^3 • We take fj (t) to be a quantity 
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FIG. 2: Distribution of events of size s, P{s). P(s) does not 
exhibit power law behavior for a = 2.5 and R = 1, but does 
so for larger R with & ~ 2. For a = 0, P{s) ~ s"*" with fe ~ 1.5 
for R > 100. Note that s can exceed A'' because a block can 
fail multiple times during an event. 

associated with block j and define 

7At)^-Jjjit')dt' (4) 
1 ^ 

If the system is ergodic, 51/ (i) oc 1/t jljj]. The met- 
ric studied in the CA models is the stress metric. In 
Fig. Owe show the inverse stress metric flj^{t) for the 
Burridge-KnopofF model as a function of the number of 
times the substrate is moved (loading times) . For R ~ 1 
the system appears to be ergodic, unlike the short-range 
CA models [ig- For a = 2.5, the system ceases to be er- 
godic for larger values of R. In contrast, for a — 0.5 the 
system remains ergodic, and the slope of fij ^ becomes 
larger as R increases. In particular, for a < 0.5 and large 
R, the behavior of f2/ is similar to the stress metric in 
the long-range CA models [l^ . 

To understand this behavior, we plot the time series 
of the mean stress </> per block just after an event 
for R = 500 (see Fig. g}. For a = 2.5 <f{t)> is quasi- 
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FIG. 3: The inverse stress metric versus the number of loading 
times. Note the different vertical scales. The system appears 
ergodic for a = 2.5 and i? = 1, but is nonergodic for larger 
values of R. For a = 0.5, the slope of Qj^{t) becomes larger 
as R increases. For a — 0.5 and R = 100, the system exhibits 
punctuated ergodicity similar to the behavior of the long- 
range CA model. 
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FIG. 4: The mean stress </> as a function of the number of 
loading times for R — 500. Note the fluctuating behavior for 
a = 0.5 and the quasiperiodic behavior for a = 2.5. 

periodic and ranges from « —0.2 to almost 1.0 with a 
mean of « 0.4. That is, small earthquakes accumulate 
stress locally and characteristic earthquakes release the 
stress globally and quasi-periodically. The periodicity of 
characteristic earthquakes for R — 1 and a = 2.5 was 
also analyzed in Ref. (3] • The distribution of characteris- 
tic earthquakes and other results 20] show that the pe- 
riodicity becomes better defined as R is increased. The 
quasi-periodic behavior in Fig. 0] is reminiscent of the 
stress versus time curves observed in laboratory experi- 
ments with rocks In contrast, for a = 0.5, <f{t)> 
fluctuates about 0.99 (see Fig.^J. The latter behavior 
is observed in long-range CA models [T^ . 

Other statistical measures of the long-range Burridge- 
Knopoff model approach those of the long-range CA 
model. For example, the mean displacement (slip) of 
the blocks during an event, <u>, as a function of s, the 
number of blocks in an event, is shown in Fig. [S] For 
large R the mean displacement is independent of the size 
of an event up to the largest events, as in the long-range 
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FIG. 5: The mean displacement <u> (slip) of the blocks as 
a function of the number of blocks s in an event. (Each block 
is counted once even if it fails multiple times.) 

CA models. The range of events over which the mean 
slip is a constant increases with R. 

In summary, we have found two possible scaling 
regimes, one of which (a — » and R » 1), is asso- 
ciated with an equilibrium spinodal critical point with 
an exponent b = 3 /2 c onsistent with the long-range CA 
models [ll|ll|ll|llli3. The power law scahng of the 
moment for a = 2.5 and R = 1 with & w 2 has been asso- 
ciated with critical behavior |^ . This latter behavior 
and the ergodicity of the stress metric for this case sug- 
gests the presence of another critical point for short-range 
systems, but this interpretation needs more investigation. 
That is, the Burridge-Knopoff model appears to have two 
scaling regimes with qualitatively different behavior. 

The apparent dependence of 6 on a suggests that larger 



system sizes as well as longer run times should be inves- 
tigated. However, because real faults are finite and the 
number of events observed is small, the a-dependence 
of b seen in the long-range Burridge-Knopoff model may 
accurately reflect the behavior of real faults. 

For a — 2.5 and large R, our results resemble those ob- 
served in laboratory experiments on rocks. As the range 
R is increased and a is made smaller, our results more 
closely resemble the long-range CA models. This wide 
range of behavior indicates that the physics of several 
models of earthquake faults ^, ^ ^, can be obtained 
from the generalized Burridge-Knopoff model with the 
appropriate choice of R and a. 

Our results imply that real earthquake faults and lab- 
oratory rocks can have different statistical distributions 
of events and different physical characteristics due to the 
details of the friction force as well as range of the stress 
transfer. This dependence means that we must develop 
ways of determining the friction force with considerably 
more accuracy if we want to understand the relation be- 
tween physical processes and observed earthquake phe- 
nomena, and if we are to be able to predict earthquakes 
or even determine if they are predictable. 
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